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SmJ • Abstract 

CO ■ An integrable discretization of the inhomogeneous Ablowitz-Ladik model with a linear 

force is introduced. Conditions on parameters of the discretization which are necessary for 

reproducing Bloch oscillations are obtained. In particular, it is shown that the step of the 

^_^ discretization must be comensurable with the period of oscillations imposed by the inhomo- 

I ' geneous force. By proper choice of the step of the discretization the period of oscillations of a 

>-i. , soliton in the discrete model can be made equal to an integer number of periods of oscillations 

^--s ' in the underline continuous-time lattice. 

o 

Q^ ■ 1 Introduction 

^' . 

During the last few years a great deal of attention has been paid to integrable discretizations of 
nonlinear evolution equations. The interest is naturally justified by needs of the computational 
physics H, 0]. One of the purposes of integrable discretization is the construction of a discrete 
analogue of a continuum model which preserves main features of the last one. This point becomes 
O ' especially important when one deals with inhomogeneous models. In that case even "the first 

,, I step" of the discretization of a one-dimensional nonlinear evolution equation, i.e. discrtization with 

respect to the spatial coordinate, may introduce qualitatively new features into the dynamics. So, 
for instance, in the case of the inhomogeneous nonlinear Schrodinger equation a constant force, 
which linearly depends on the spatial coordinate, results only in the renormalization of the phase 
cH I and velocity of the one-soliton solution [^ while the same force leads to oscillations of solitons of 

the inhomogeneous Ablowitz-Ladik (AL) 0] model: 



^ 



iqn + (1 - gn7-„)(g„-i + g„+i) + 2xnqn == (1) 

- ir-n + (1 - (?„r„)(r„_i + r„+i) + 2xnr„ = (2) 

(here r„ = ±g„, x is a real constant, which from the physical point of view determines the strength 
of the linear force, a dot stands for the derivative with respect to time, and a bar stands for the 
complex conjugation) |5| ||, 0]. Periodic dependence on time is a property of any solution of (|l]), 
(0) and it is caused by the discreteness. In the case of the one-soliton solution which reads 



_sinh(2w)^^__^.[<I.(t)-2nx(t-to)] 

cosh[2^ 
where 



^" " cosh[2nw-X{t)-Xo] ^' 



$(i) = _ cosh(2w) sin[2x(i - to)] (4) 

X 



X(t) = - sinh(2u;) cos[2x(i - ^o)] (5) 

X 

w, io, and Xq are real constants, the soliton dynamics has deep analogy with the well known Bloch 
oscillations of an electron in a lattice potential affected by a constant electric field (due to this 
reason such behaviour is referred to as Bloch oscillations g] ) . As it follows from (^ the period of 
oscillations is given by tq = tt/x- 

The phenomenon of Bloch oscillations becomes especially interesting if one looks for the pos- 
sibility of integrable discretization of the inhomogeneous AL model with respect to time. Indeed, 
Bloch oscillations are characterised by the additional temporal scale, tq. This scale is determined 
by the strength of the force and must lead to some constrains on the step of discretization. Thus 
the purpose of the present communication is to introduce integrable discretization of the model 
(0) , (0) and to obtain conditions on parameters of the discretization which preserve the effect of 
Bloch oscillations in the discrete scheme. 

2 Integrable discretization 

At X = system (P, (|]) transforms to the conventional AL model which discretization is well 
known |g, g|. In particular, it can be achieved by using the discrete analogue of the zero-curvature 
condition 

U{n, t + h)V{n, t) = V{n + 1, t)U{n, t). (6) 

In the case x ^ the same condition involving [/-matrix as follows 

and ^-matrix having the elements 

Vii=i~- hao + hai (^A2e-*x('*+"'') - A(n, t)) + h (^ Je^^^'*-"'') - 52q{n, t + h)r{n - 1, i)) A(n, t) 

yi2 = h ('aiAe-*>^(*+"'')g(n, t) - ^e'^(*+("-i)'')g(n -l,t + h)\ + 
h (52Ae-^'^(*-("+i)'')g(n, t + h) - ^e'^'-'~''''\{n - 1, i)) A(n, t) 

V^2i = h ('aiAe-*x(*+("-i)''V(n -l,t + h)- y e*x(*+"''V(n, i)) + 
h ('^2Ae-*x(*-"'^V(n - 1, t) - ^e*>^(*-("+i)''V(n, t + h)] A{n, t) 

V22=i + hSo - hSi Ux(^t+"h)/x^ - D{n,t)) - h ((JiA^g-'^^^t-n/^) _ q,^^^^ _ l,t)r{n,t + h)) A{n,t) 

where A is a spectral parameter, aj and dj are parameters and ft. (ft. > 0) is a step of the discretiza- 
tion, results in the system 

ih'^[q{n, t + h)- q{n, t)] = 5iq{n -l,t + ft)e*'''^("-2) _ ^^^(^^ t + h) ~ 

aoqin, t) + aiq{n + 1, <)e-'''x("+i) - 5iq{n, t + h)D(ji, t) - aiq{n, t)A{n + 1, t) 

a2q{n ~l,t) [e-'''x("+i) - q{n, t + h)r{n, t + h)]A{n, t) + 

+52q{ii + 1, t + ft)[e^''x("+2) - q{n, t)r{n, t)]A{n + l,t), (8) 



~ ih-^[r{n,t + h) - r{n,t)] ^ air{n - 1, t + /i)e"*'''^(""2) -~aor{n,t + h) - 

Sor{n, t) + 5ir{n + 1, t)e*'^'^("+i) - air{n, t + h)A{n, t) - (5ir(n, t)D{n + 1, t) + 

52r{n - 1, t)[e'''x("+i) - q{n, t + /i)r(n, t + /?,)]A(n, t) + 

a2r{n + l,t + h) [e-*''x("+2) _ ^^-^^ ^y^^^ ^^^^^ ^^^^-^ 



(9) 



Qfl 



ai [A(n + l,f)-A(n,i)e-*''x] - /i^i(i - /lao) (l - e^^'^x) = 

r-(n, t)q{n + 1, t)e-^''x("+i) - r(n - 1, f + h)q{n, t + h)e-'''^^''-^A + 

52r{n - 1, t)q(n, i + h)K{n, t) (e"*'''^ - e*''"^) - 

(52<?(n + 1, i + /i)r(n, t) (l - e'^^'''-'+^A A(n + 1, t) 



(10) 



(5i [D{n + l,t)-D{n,t)e'^^] + /i-i(i + Mq) (l - e"*^) = 
(Si [r(n, t + /i)g(n - 1, i + /i)e^''x(»-i) - r{n + 1, i)g(n, t)e^''x(»+i)l + 
a2r(n, i + h)q{n - 1, t)A(n, i) (e*'*'^ - g"*'*'^") + 
a2r(n + 1, t + ;i)g(n, t)A(n + 1, t) ('e-'''x("+2) - l) (11) 

A(?i, t) [1 - g(n, t + h)r{n, t + h)] = A(n + 1, f) [1 - q{n, t)r{n, t)] (12) 

Then the discrete analogue of the AL model (|^), (g) is obtained from (g)-(|lj) by means of the 
reduction 

r{n,t) ^ ±q{n,t) (13) 

which requires the foUowing relation among the parameters: aj = dj. 

In order to define solutions of (p|)-(p^ one has to fix boundary conditions for r{n,t), q{n,t), 
A{n, t), and D{n, t). In what follows we deal only with the case of zero boundary conditions 



lim q{n,t) ~ lim r{n,t) ~ 

n — ^itoo n — >iboo 



(14) 



which allow existence of "bright" solitons. Hence it will be assumed that r{n,t) = -~q{n,t). 
Respectively we have to require 



lim A{n,t) = Ao{n,t) 

\n\ — >oc 



-i + hag 
hai 



(^ 



-ihxn 



1 



lim D{n, t) = Do{n, t) = '-±^ {^^^^ - l) . 



(15) 



(16) 



Notice that dl3) , dl9) transform to the zero boundary conditions in the case of the homogeneous 
AL model (x = 0)^. 

Eqs. (p^), (O) subject to dlq), (llq) allow one to express Ain^t) and D{n,t) through q{n,t) 
and r{n,t): 



A{n, t) = Ao(n, t)+a^^Y^ fA{n - k, t)e 



— ixhk 



fc=l 



Din,t)^Doin,t) + 5^'J2fD{n-k,t)e'^'''' 



(17) 



(18) 



k=l 



Here /yi(n,i) and foin^t) stand for the right hand sides of the equations ( |lO| ) and (|T|), corre- 
spondingly. It is to emphasised that formulae ( p7| ) and ( p^ ) do not give yet explicit solutions for 



A{n,t) and D{n,t) and after substitution into (p[), (0) represent a source of nonlocality of the 
discrete scheme. 

As it has been shown in l6l 0] a convenient approach to treat inhomogeneous discrete models 
is the use of the gauge transformation which allows one to restrict the study only to the temporal 
behaviour of the scattering data. The gauge transformation in the discrete case is given by 

tj{n,t) ^ G{n + l,t)U{n,t)G-\n,t) (19) 

V{n,t) =G{n,t + h)V{n,t)G-^{n,t) (20) 

By choosing G{n,t) — explixnta^}, where CT3 is the Pauli matrix one reduces U{n,t) to the 
form which corresponds to the f/-matrix of the underline homogeneous model (i.e. to the form 
which does not have explicit dependence on time and can be obtained from m) by the replacement 
exp(zxt) K^ 1). Then the dependence of the transfer matrix T(i), associated with U{Ti,t), on the 
discrete time is governed by the equation 

T{t + h) = VhT{t)V,:' (21) 

where Vh is a diagonal matrix, Vh =diag(6'i(A, t), 6'2(A, i)), with the elements 

ei{X,t) =i-haa + /iA^aie-2«xt + hX-^a2e^'^K (22) 

6'2(A, t)^i + hSo- /lA^fJae-^^xt _ hX-^Sie^'^K (23) 

Let us now assume that t — mh and tti = 0, 1, ... (it is straightforward to generalise the results 
to the case t — mh + Iq where to is an arbitrary real constant playing the role of initial moment 
of time). Then the element T^^^^ of the matrix T{t) does not depend on m (or t) while for 
r(i2)(i) = 5^ one obtains 

n=Q 

where 

^i„,{X) = 0i{X,mh)/02{X,mh). (25) 

In the case of solitonic solutions (p4) formally solves the discrete Cauchy problem since it 
defines dependence of the scattering data on time and the solution of the eigenvalue problem for 
the matrix U{n,0) is well known B. Below we concentrate on some "physical" consequences of 
that result. 

3 Bloch oscillations in the discrete model 

As it has been mentioned in the Introduction each step of discretization can introduce new features 
into the dynamics (even in cases of integrable models). One of such features is the oscillatory 
behaviour of the solutions of the AL model affected by the linear force (Bloch oscillations). Now we 
address to the question whether it is possible to preserve such evolution subject to the discretization 
with respect to time. 

To this end we take into account that periodic behaviour means that there exists a positive 
integer M such that 

m+M 

n Mn(A) = l (26) 

n— 771+1 

for any A (which can be considered, say, inside the unit circle on the complex plane) and any m. 
The period r of oscillations is then given by r = Mh (evidently M is considered to be the smallest 
possible integer). 



The discretzation of the homogeneous AL model is a three parametic one (this, in particular, 
allows one to represent it in a form of a product of the local maps [H). In the inhomogeneous case 
the imposed conditions lead to constrains on the parameters. To find them we first consider the 



limit |A| -> (or |A| -^ oo). Then from (|2^), (g3|), (^, and (|26|) one finds that there must exist 
relations 

ai = -a2 0f^ + 4"=2^^, Z = 0,l,...,Af-l (27) 

where 01^2 and 4>i^2 are real prameters connected to ai.2, ai^2 = Q1.2 exp(0i_2), and the upper 
index has been attributed to the "quantized" phases. Next, the independence of (Pq ) on m implies 
Mm — f^m+M which means that 

xhiM = Til (28) 

where I is a positive integer. In other words the step of the discretization is not arbitrary (the 
subindex I is introduced to label different discrete values of h) . Physical sense of the last require- 
ment is quite transparent. Recalling that the period of Bloch oscillations in the continuous-time 
model is given by tq = tt/x, one concludes that (|2^) means that the period of the Bloch oscilla- 
tions in the discretized model, r, is I times bigger than the period tq: r = Itq and the number / 
is related to the chosen step of discretization h. On the other hand rewriting ( P8[) as hi = Itq/M 
one can interpret it as a condition for the discretization step to be comensurable with the period 
of oscillations. As it is evident, for the direct coincidence of the result obtained on the discrete 
lattice with its continuum counterpart one must let Z = 1. Below we concentrate on this case. 
Then /in (A) takes the form 

1 + ap*(^'+'''-")A2 - ap*(^'-T''")A-2 



where 



r ^ A ^ (A'^) J^A "^^^ 



a — \hai/{i — hao)\ and ^o = arg(i — hao). 

Now we consider the unit circle where A^ = exp(i7/;) (V' being real). Then one can find two 
possibilities to satisfy the requirement (|26|). One simplest solution corresponds to M even and 
(j)o — ttI/M + tt/2 + Tip {p is an integer). Then /i„(A) = exp{27ri(Z/M — 1/2)}. By the direct 
algebra one ensures that this is the degenerated case, when the limiting transition /i — > results 
in a trivial linear equation instead of the AL model. A nontrivial and physically relevant solution 
corresponds to the case when M ~ 4N {N is an integer) and (po — (pi.p = ttI/M + irp (in that 
case iinfJ-n+NtJ'n+2NfJ'n+3N — !)• Then /iri(A) which determins evolution of the one-soliton solution 
associated with the eigenvalue Ai = exp(— w -I- i6) is given by 

^"(^^) = ' l + 2(-l)Pasinh[2^-z(7.„ + 2^)] °"P (.'"^mJ ^^^^ 

Let us illustrate the discrete-time dynamics on example of the one-soliton solution. We assume 
that ( p7| ) holds. For the sake of simplicity we let 0i = —(j)2 = 7r/2 and 0o = 0. Then the one-soliton 
solution of (g), (g) can be written down in the form (recall t = mh) 

q^^\n,m) = -&\n,m) = ^inh(2^) ^,,(„,„, 

cosh(2nw — Xm) 



where 

m — 1 



4>(n,TO) = \^ arctan 



fc=0 



4a cos {2xkh - 20) cosh(2w) 



1 -I- 2a^ cos {Axkh - A9) + 4a'^ cosh(4w) 



2ixnmh (32) 



fe=0 



2c? (cosh(4w) + cos (4xfcft. - 46*)) + 1 - 4a sinh(2u;) sin {2xkh ~ 29) 



2a2 (cosh(4w) + cos (4xfc/i - 46*)) + 1 + 4asinh(2u;) sin (2xA:/i - 20) 



(33) 



Comparing these formulae with (|j), (||) one can see that choosing a = h/2 the last ones are 
obtained by the limiting transition /i — j 0. 

To be more specific we concentrate on X^ which describes evolution of the centre of the soliton. 
To this end we represent h = ttS,/{Mx) (the so chosen step of the discretization satisfies (GS) when 
^ is integer) . The results are summarized in Fig. 1 for three values of the parameter ^ displaying 
different situations: (a) When ^ = 1 the discrete model exactly reproduces Bloch oscillations of the 
continuous-time model, (b) At ^ = a/3 the evolution of the discrete model is not periodic (notice 
that lines in the figure are used for the sake of convenience of presentation: the truth trajectories 
are sets of points). However, by considering analytic continuation of the solution the periodicity 
can be considered between discrete time steps M. (c) At ^ = 2 the period of soliton oscillations 
is two times more than the period of Bloch oscillation of the AL model which is obtained by the 
limiting transition /i — > 0, M — > oo with hM — tt/x (it is to be mentioned that the minima of the 
curve (c) about h — 48 and h — 145 are not numerical zeros). 

Author acknowledges support from FEDER and Program PRAXIS XXI, grant N° PRAXIS/2/2. l/FIS/176/94. 
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Figure caption 

Trajectory of the soliton centre corresponding to different steps of discrete time (a) h = 1/97, 
(b) h — \/3/97, and (c) h = 2/97. Other parameters are as follows: M — 97, w — 0.5, x = ^r, 
a = 0.1. 
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